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The investigation of the dynamics of quantum many-body systems is a concerted 
effort involving computational studies of mathematical models and experimental 
studies of material samples. Some commonalities of the two tracks of investiga- 
tion are discussed in the context of the quantum spin dynamics of low-dimensional 
magnetic systems, in particular spin chains. The study of quantum fluctuations in 
such systems at equilibrium amounts to exploring the spectrum of collective exci- 
tations and the rate at which they are excited from the ground state by dynamical 
variables of interest. The exact results obtained via Bethe ansatz or algebraic 
analysis (quantum groups) for a select class of completely integrable models can 
be used as benchmarks for numerical studies of nonintegrable models, for which 
computational access to the spectrum of collective excitations is limited. 

1 Introduction 

Notwithstanding the fact that experimental and theoretical studies of con- 
densed matter systems are fundamentally complementary to each other, they 
share important features, which we wish to illuminate here in the context of 
quantum many-body dynamics. The object of study is a sample for experi- 
mental probes and a model for computational probes. Whether the model is 
regarded as a mathematical idealization of a real chunk of matter or whether 
the sample is viewed as a physical realization of a system defined with math- 
ematical rigor may be a matter of perspective, but the influence exerted by 
theory and experiment is mutual. Any observation of note calls for an ex- 
planation, while any prediction of substance brings into motion attempts at 
verification. 

When performing experimental or computational probes in quantum 
many-body dynamics, the goal of researchers is very much the same, namely 
trying to make sense of what, in general, is a jumble of fluctuations to which 
the fundamental degrees of freedom contribute in various combinations and 
configurations. Experimental probes begin with a measurement and compu- 
tational probes with a calculation. The results of either probe then require an 
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interpretation in terms and concepts that are common to both approaches. 

Qualitatively, the specification of a quantum many-body system involves 
information on composition, interaction, and environment. Theoretically 
this information is encoded in the Hamiltonian and in the density operator, 
whereas experimentally part of this information is manifest or hidden in the 
sample and other parts are controlled by the instrumental setup. 

Any system thus specified is subject to fluctuations, of which we distin- 
guish three kinds. Depending on the circumstances, each kind may govern 
the dynamical properties of the system. Thermal fluctuations are likely to be 
dominant in any statistical ensemble at elevated temperatures. Parametric 
fluctuations are manifestations of quenched random inhomogencities in com- 
position (e.g. dilution) or interaction (e.g. random bonds or fields), leaving 
distinct marks in dynamical properties, especially at low temperatures. 

In the absence of thermal and parametric fluctuations, quantum fluctu- 
ations remain present. They are a direct consequence of the (autonomous) 
quantum time evolution and the time-delayed projections that are part of any 
dynamical probe, be it experimental or computational. No zero point motion 
exists in classical Hamiltonian systems. Here the ground state corresponds to 
a stable fixed point in phase space, and all dynamical variables are constant. 

Consider a (non-random) quantum many-body system in a pure quan- 
tum state. The quantum fluctuations then depend on three quantities: (i) 
the Hamiltonian H, (ii) the state, which we take here to be the ground state 
\G), and (iii) some dynamical variable, here expressed by operator A. The 
Hamiltonian governs the dynamics deterministically for as long as the sys- 
tem can be regarded in isolation. It makes no difference whether this time 
evolution is viewed in the Schrodinger or Heisenberg representation: 



ih-m=HW) W))=e- iHt/ », \^)=A\G) (1) 



The statistical element is introduced when we subject the time evolved state 
\ip(t)} or, equivalently, the time evolved dynamical variable A(t), to an ob- 
servation. In practice, this step involves the observation or evaluation of a 
dynamic correlation function, which can then be viewed either as the pro- 
jection of the time evolved state \ip(t)) or as the product of the dynamical 
variable A(t) at two different times evaluated in the stationary state \G): 



[A,H] 



A(t) = e 



,iHt/hj^ e -iHt/h 



(2) 




(3) 
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The sum total of quantum fluctuations in a typical many-body system 
contains many different dynamical modes no matter in which state the system 
happens to be. Only once we pick a set of dynamical variables (one variable at 
a time) do we begin to gain insight into the role of specific modes taken from a 
multitude of quantum fluctuations. Criteria for choosing dynamical variables 
include experimental accessibility, computational amenability, or any specific 
purpose like the study of ordering tendencies or phase transitions. 

On the microscopic level, dynamical probes are bound to be somewhat 
heavy-handed. No observation without perturbation! However, in all situa- 
tions considered here, we shall assume that the interaction hetween the probe 
and the sample takes place in the regime of linear response.a This assumption 
is quite realistic for neutron scattering under most circumstances. Consider 
the weakly perturbed Hamiltonian H(t) = Hq — b(t)B, where the external 
field b(t) couples to the dynamical variable B of the system Hq. The linear 
response of any other dynamical variable A to that perturbation, 

/+oo 
dt'xAB(t-t')b(t% (4) 
-co 

is then determined by Kubo's formula for the response function, 

X AB (t-t') = ^Q(t-t')([A(t),B(t')}) . (5) 

Its Fourier transform, the generalized susceptibility xab{u + has (for 
e — > 0) a symmetric real part and an antisymmetric imaginary part: xab (w) = 
x'ab^) + ^Xab^)- The latter is related to the Fourier transform of the 
correlation function (A(t)A), named the structure function 

SaaM = [ + °° dte^(A(t)A) = f*^ , (6) 

via the fluctuation-dissipation theorem. The spectral representation of this 
function, in particular the expression resulting in the limit T — > 0, where 
thermal fluctuations are absent, 

S aa (lo) 2ttJ2 \(G\A\X)\H (w - (E X - E G )/K) , (7) 

A 

demonstrates that observing quantum fluctuations in the linear response 
regime means observing collective modes and their transition rates from the 
ground state. In the following, we look at quantum fluctuations and the asso- 
ciated collective modes from this particular angle for a number of situations 
as investigated by a variety of methods]^ 

"From here on we set h = 1 to simplify the notation. 
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2 Magnons Excited from Valence-Bond- Solid State 

Consider the one-dimensional (ID) s = 1 Hciscnberg model with bilinear and 
biquadratic exchange 

N 



At T = it has a phase with ferromagnetic long-range order (LRO) (— tt < 
9<-3tt/4 or tt/2<9<tt), a phase with dimer LRO (-3tt/4 < 9 < -tt/4), the 
Haldane phase with hidden topological LRO (— 7r/4< < 7r/4), and an obscure 
phase (tt/4 < 8 < tt/2) that was named trimerized.0 These ordering tendencies 
help us choose suitable dynamical variables when we explore the spectrum of 
collective excitations. All dynamical variables used for this purpose will have 
the form of fluctuation operators, 



where the operator A n acts locally at or near lattice site n. The associated 
structure function (ub is called the dynamic structure factor. 



where the sum runs over the dynamically relevant excitations |A) with en- 
ergy uj\ = E\ — Eg and spectral weig ht W£ = 2tt|(G|F^|A)| 2 . To calculate 
Saa(Q,w) for Hg we can employ special methods at the parameter values 
9 = ±7r/4 where it is completely integrableS or general methods otherwise. 
One suitable general method in the present context is the recursion method 
in combination with a finite-size continued-fraction analysis.Q 

Here we introduce four different fluctuation operators Fq. The local 
operators A n from which they are constructed are listed in Table [l| 

• The spin fluctuations, probed by Fq, represent Neel order parameter fluc- 
tuations for Q = tt. They are expected to be strongest at the critical point 
6 = — tt/4, where the Q = tt excitations are gapless. 

• The dimer fluctuations, probed by Fq, are also expected to be strongest 
(for Q — tt) at 6 = — tt/4, which marks the onset of dimer LRO. 

• The trimer fluctuations, probed by Fq, are constructed from projection 
operators P£ onto local trimer states |[1,2,3]). The state |[1,2,3]) happens 
to be the ground state of Hg with N = 3 for arctan | < 9 < tt/2, which was 
interpreted as suggesting that a trimerized phase might exist for N — > oo. 
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(10) 



• The center fluctuations, probed by Fq, are constructed from a modified 
spin operator and tune into existing period-three (+,0,—) or (—,0,+) pat- 
terns of local spin states. Finite- N data indicate that for 7r/4 < 8 < tt/2, the 
fluctuations probed by F^ w , 3 are the strongest of all the ones listed. 



Table 1. Local operators A„ which define the four fluctuation operators Fq via (g) and 
the four order parameters Pa via (Jill) . 
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The order parameters associated with the four fluctuation operators de- 
fined by Eq. (^) and the entries of Table [l] can be written in the form 



p A = ^Y. elQAnA ^ ( n ) 

n=l 

where Qs = Qd = 7r, and Qt = = 27r/3. Each order parameter Pa has a 
set of eigenvectors k = 1, 2, . . . , -K" which represent the associated LRO 

in its purest form. The degree of degeneracy is K = 2 for the Neel and dimer 
states, K = 3 for the trimer states, and K = 6 for center states.B 

The physical vacuum chosen here is very unlike any of the states 
Within the Haldane phase, at the parameter value Ovbs = arctan^, the 
ground state of Hg is a realization of the valence-bond solid (VBS) wave 
functionjj which is non-degenerate and in which the Neel, dimer, trimer, and 
center ordering tendencies are all imperceptibly weak. In the VBS state, the 
spin 1 at each lattice site is expressed as a spin-1/2 pair in a triplet state. The 
singlet-pair forming valence bond involves one fictitious spin 1/2 from each of 
two neighboring lattice sites. The VBS state, which has total spin St = 0, 
can then be regarded as a chain of valence bonds linking successive spin-1/2 
pairs in this manner. 

The topological LRO present in the Haldane phase and known to be 
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strongest in the VBS state provides an environment, where a specific kind of 
elementary excitations can propagate freely and where the associated station- 
ary states form a branch with well defined dispersion. We now probe these 
elementary excitations and composites thereof from different angles by the 
four fluctuation operators Fq defined in Table [l]. 

Only two of the Fq are fully rotationally invariant, [Fq, Sj] — [Fq, Sj] =0 
for i = x,y,z. This produces the selection rules ASt — in the dynamic dimer 
and trimer structure factors. The corresponding selection rules for the dy- 
namic spin and center structure factors are ASt = 0,1 and ASt = 0,1,2, 
respectively, with the further restriction that transitions between singlets 
(St = states) are prohibited. At the VBS point, Sdd{Q,u) and Stt(Q,u) 
thus couple exclusively to the St — excitation spectrum, and SssiQj^) 
exclusively to the St = 1 excitation spectrum, whereas Szz(Q,<+>) couples to 
the St = 1 and St — 2 spectra. 

In Fig. [I] we display lo^ versus Q of the dynamically relevant spin, center, 
dimer, and trimeiiexcitation spectra as obtained from the finite-size continued- 
fraction analysis The filtered access to the spectrum afforded by the four fluc- 
tuation operators gives us valuable clues about the nature of the elementary 
excitations that thrive in the VBS environment. The low-frequency region at 
Q > 7r/2 in the spin and center spectra is dominated by a branch of states 
with St = 1. They carry more than 95% of the spectral weight in Sss{Q,u) 
and Szz(Q,u)- These triplets, which have been named magnons, remain in- 
visible in the dimer and trimer spectra. Only composites of the magnon states 
may be observable in Sdd(Q,u) and Stt(Q,^)- 

A dispersion of the general form com (Q) = J(a + b cos Q) for the 1-magnon 
branch can be inferred from the single-mode approximation of Sss(Q, co)nu 
Under the assumption that in the VBS environment the magnons are weakly 
interacting point particles, we can expect the existence of three kinds of 2- 
magnon scattering states formed by pairs of 1-magnon triplets: states with 
St = 1, which contribute to Sss(Q,u) and Szz(Q,u), states with St = 0, 
which contribute to Sdd(Q,^) and Stt{Q,oj), and states with St=2, which 
are observable in Szz{Q, w) only. Free 2-magnon states form a two-parameter 
continuum L02M(k,Q) = u>m(Q/2 — k) +wa/(Q/2 + k) in (Q, w)-space. The 
resulting continuum boundaries are u>±(Q) = 2J[a ± 6cos(Q/2)]. 

The predicted coalescence of the 2-magnon continuum into one spectral 
line at Q = ir follows from the symmetry property uim{Q)+^m(^—Q) = const 
of the magnon dispersion. Interestingly, the finite-TV dimer spectrum does 
indeed collapse into a single spectral line at the iV-independent excitation 
energy ujd = \/T0J, which carries all the spectral weight in Sdd(Q,u). 

We now use the exact 2-magnon excitation energy, lu± (jr) = lob , and the 
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Figure 1. Frequency u>jr versus wave number Q of the finite-iV excitations [N = 12 (*), 
N = 15 (O), iV = 18 (B)] which carry most of the spectral weight in the T = dynamic 
structure factors SaaCQi^) f° r the spin, center, dimer, and trimer fluctuations of Hg with 
J = 1 at 8 = arctan 1 . The three sizes of symbols used distinguish excitations with relative 

spectral weight w£ = W^/(F^F^) in the ranges w£ > 0.5 (large), 0.5 > w£ > 0.1 

(medium), 0.1 > w£ > 0.001 (small). The dashed lines represent 1-magnon dispersions 
and the solid lines 2-magnon continuum boundaries as explained in the text. 

extrapolated value, ojm{^) = 0.66433(2) J, of the 1-magnon excitation gap to 
fit the parameters a, b. The resulting values, a ~ 1.581,6 cp. 0.917, used in 
wm(Q) and ijJ±{Q) yield the dashed and solid lines in Fig. pl.El 

Both the 1-magnon dispersion and the 2-magnon spectral threshold are 
in very good agreement with all finite-A^ data shown. Hence the magnon 
interaction is very weak at the bottom of the 2-magnon region in all three 
St subspaces. The finite- N data spilling out on top of the shaded areas in 
Fig. ^ suggest that at higher energies, the magnon interaction is repulsive, 
more strongly so in the St = 0, 2 subspaces than in the St = 1 subspace. 

This raises the possibility that bound 2-magnon states split off the top of 
the 2-magnon continuum of 2-magnon scattering states in the St = 0, 2 sub- 
spaces. The comparison of panels (a) and (b) at frequencies 3 < u/J < 5 in- 
deed suggests that the dynamically relevant finite- TV excitations are arranged 
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in contrasting patterns. In panel (a) we have an arrangement of points which is 
typical of a two-parameter continuum. As N increases, more points are added 
and spread roughly evenly along the frequency axis. In panel (b), by con- 
trast, the data points are arranged in branches with an almost TV-independent 
separation, which is characteristic for branches of bound states. The spin-1 
compounds Ni(C 2 H$N 2 ) 2 N0 2 Cl0 4 (NENP) and Ni(C 3 H 10 N 2 ) 2 N 3 {ClO 4 ) 
(NINAZ), while not physical realizations of Hg directly, nevertheless realize 
situations where the spectrum of collective excitations as probed by inelastic 
neutron scattering^ shares major features with the magnons excited from the 
VBS state. 



3 Magnons Excited from Ferromagnetic State 

We now turn to a more quantitative discussion of the difference between scat- 
tering states and bound states in the context of completely integrable situa- 
tions, namely for the ID s — | Heisenberg ferromagnet: 

JV 

Hp = — J S„ • S n+ i. (12) 

The ground state is (iV + l)-fold degenerate. We select one of the ground- 
state eigenvectors, \F) = | ft " " T)i as the physical vacuum for an explo- 
ration of collective excitations. As in the VBS case discussed in Sec. ||, the 
spin fluctuation operator probes transitions between the vacuum state and a 
branch of 1-magnon states. However, here the magnon dispersion is gapless: 
E-E F = J(l-cosk). 

Unlike in the VBS case, the 1-magnon states of Hp are located in a 
separate invariant subspace. Only the 1-magnon states can contribute to the 
spin fluctuations. In the expression 

1 N 

\k) = S^\F), S- = _£y*»S- ( 13 ) 

v n— 1 

for the 1-magnon eigenvectors, the spin fluctuation operator plays the role of 
a magnon creation operator. That is not to say is a true magnon creation 
operator. Multi-magnon superpositions, i.e. the states S% ■ ■ ■ \F), are a 
redundant set of non-orthogonal and non-stationary states, which are subject 
to two kinds of interactions]!! 

• The kinematical interaction is caused by the restriction on the number of 
reversed spins at one lattice site. 
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• The dynamical interaction is caused by the off-diagonal part of Hp in the 
basis of multi-magnon superpositions. _ 
This distinction is quite natural in the framework of the Bethe ansatza 
as we shall see. The Bethe ansatz is an exact method for the calculation 
of eigenvectors and eigenvalues of completely integrable quantum many-body 
systems. The Bethe wave function of any eigenstate in the subspace with 
r = N/2 — S?p reversed spins relative to the magnon vacuum, 

IV;} = a(n 1 ,...,n r )S- 1 ---S- r \F). (14) 

l<ni<...<ri,-<7V 

has coefficients of the form 

( r i r \ 

0(711,... ,n r ) = exp ij^^rij + -J^^Vj (15) 

ves r \ 3=1 i<j J 

determined by r magnon momenta ki and one phase angle <?y = — 6ji for 
each magnon pair. The sum V € S r is over the permutations of the labels 
{1, 2, . . . , r}. The consistency requirements for the coefficients a(ni, . . . , n r ) 
inferred from the eigenvalue equation H\^p) = E\t/j) and the requirements 
imposed on the same coefficients by translational invariance can be cast in a 
set of equations for the momenta ki and phase angles Oij\ 

9- k k 

2cot^- = coty - cot—, ^=27^ + 2^%. (16) 

Every solution of these equations is specified by a set of r Bethe quantum 
numbers Xi G {1,2,... , N — 1}. Given a solution, the energy and wave 
number of the state it describes are 

r 2tt r 

E-E F = J^2(l-cosk j ), k = —J2 X *- ( 17 ) 

j=i i=i 

In the subspace with r = 1, we thus recover all N 1-magnon states, one for 
each of the allowed values of Ai. There exist N(N +l)/2 distinct 2-magnon 
superpositions of the 1-magnon states thus identified. However, this set of 
states must be accommodated in the r — 2 subspace, whose dimensionality 
is only N(N- 1)/2. Inspection shows that of the N(N +l)/2 pairs of Bethe 
quantum numbers in the allowed range 0<Xi<\2<N~l,N pairs do 
indeed not produce a solution of ([h]) . The missing solutions are a consequence 
of the kinematical interaction between magnons. 

The consequences of the dynamical 2-magnon interaction are illustrated 
in Fig. |^, which shows the complete r = 2 spectrum (E — Ep) / J versus k 
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for k > and N = 32.0 Also shown in Fig || are the (fictitious) 2-magnon 
superpositions, where the fci,fc2 in Eq. ( |l7| ) are replaced by all combination 
of 1-magnon wave numbers. There are three classes of states. 




Figure 2. Ex- 

citation energy 
(E — Ep)/J versus 
wave number k 
of all eigenstates 
in the invariant 
subspace with 
r = 2 overturned 
spins for a system 
with N = 32. Also 
shown are the 
corresponding data 
for free 2-magnon 
superpositions. 



The class C\ contains N states for which one of the two Bethe quantum 
numbers is zero. This means that one of the magnons has zero wave number. 
Its effect is a slight rotation of the magnon vacuum, in which the other magnon 
is as free to propagate as in the original vacuum. There is no dynamical 
interaction. All states in this class are, effectively, 1-magnon states. 

The class C2 of states has nonzero Bethe quantum numbers A2 — Ai > 2. 
There are N(N—5) /2+3 such pairs. All of them yield a solution with real fci , , 
which makes them 2-magnon scattering states. A measure of the dynamical 
magnon interaction in Fig. ^|is the vertical displacement of any true scattering 
state (o) from the nearest free-magnon pair (+). As N increases, the energy 
correction diminishes for all class C2 states and vanishes in the limit N—^oo. 
The 2-magnon scattering states and the free 2-magnon states then form two- 
parameter continua with identical boundaries E — Ep = 2J\V ± cos(/c/2)]. 

The class C3 of states has nonzero Bethe quantum numbers which either 
are equal A2 = Ai, or differ by unity A2 = Ai + 1. There exist 2N — 3 such 



pairs, but only N— 3 pairs yield solutions of (16). For the class C3 states, the 
effects of the dynamical magnon interaction are much more prominent, and 
the interaction energy does not disappear when N — > 00. In Fig. |2|, these states 
form a branch of 2-magnon bound states with dispersion E—Ep = i J(l— cos k) 
below the continuum of 2-magnon scattering states. 

The bound state character of the class C3 states manifests itself in the 
enhanced probability that the two flipped spins are on neighboring sites of 
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the lattice. This property of the wave function is best captured in the weight 
distribution |a(m,n2)| of basis vectors with flipped spins at sites n\ and rii. 
In Fig. H we have plotted |a(m,pa)| versus n 2 — n% for a sequence of class C3 
states between k = and k = irzQ The distribution is peaked at 712 — n\ = 1. 
Its width is controlled by the imaginary parts of k\ , ki in ( |l5| ) . The smallest 
width is observed in the bound state at k = it. In this case, all coefficients 
with n 2 ^ n\ + 1 are zero, which implies that the two down spins are tightly 
bound together and have the largest binding energy. 

The width of the distribution |a(ni,rt 2 )| increases as k decreases, and the 
binding of the two down spins loosens. For finite N, the Bethe ansatz solutions 
switch from complex to real when the distribution has acquired a certain 
width. In scattering states the distribution |a(ni,n.2)| is always broad and 
tends to oscillate wildly. Some scattering states have a smooth distribution 
with a maximum for n 2 = n\ + N/2, when the two down spins are farthest 
apart. The formation of bound states and scattering states of elementary 
excitations exist in many different contexts. But only in rare cases such as 
this one can they be investigated on the level of detail presented here. 

4 Spinons Excited from Spin-Fluid State 

Turning our attention to the Heisenberg antiferromagnet, we consider the 
Hamiltonian Ha = —Hp with Hp defined in Eq. All the eigenvectors 



zuoz2.1: submitted to World Scientific on February 1, 2008 11 



remain the same, but the energy eigenvalues have the opposite sign. The 
magnon vacuum \F) now is at the top of the excitation spectrum. The ground 
state | A) of Ha is located in the invariant subspace with Sj, = 0. This 
subspace also contains the two Neel states \Afi) = | tit • ■ • I), I A/2) = | ITI 
• • • j), which are not eigenstates of Ha- In the framework of the Bethe ansatz, 
\A) can be obtained from \F) by exciting r = N/2 magnons with momenta fcj 
and (negative) energies —.7(1 — cos fcj). The Bethe quantum numbers for this 
state are {Aj}a = {1, 3, 5, . . . , N — 1}. 

For reasons of computational convenience we rewrite the Bethe ansatz 



equations (16) in terms of the variables Zi = cot(fe,/2): 



N arctan Zi = nli + ^ arctan ( — - j , i = 1, . . . , r. (18) 

The associated Bethe quantum numbers —N/2 < 7j < N/2 are integers for 
odd r and half integers for even r. The relation between the sets {Xi} and {Ii} 
depends on the configuration of the solution {2^} in the complex plane. Given 
the solution {zi, . . . , z r } of Eqs. ( |1S| ) for a state specified by . . . , 7 r }, its 
energy and wave number are 

E-E F ^ 2 , 2tt ^ T 

— — = - ^TT^F' * = -r-^E*. ( 19 ) 

2—1 A Z— 1 

with Ep = JN/A. For states with real {z^, Eqs. ( |l8| ) can be converted into 
a convergent iterative process: 

^ i + l£arctan[^iH) (20) 

with starting values z,- 1 "* = nl^/N. For the ground state |A) we have 

( N 1 N 3 N 11 
Wx={- T + 5 , -4+2,-. 4-5}. (2D 

High-precision solutions {zi} can be obtained with little computational effort. 
The ground-state energy per site for ./V = 4096, for example, reproduces the 
exact result (Ea — Ep)/JN = — In 2 of the infinite chain to within 1 part in 
a millionO The distribution of magnon momenta in the ground state \A) is 
broad and peaked at fcj = 7r: 



Poiki) 



> sin — cosh — cot 
2 



(22) 
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This state, which has obviously a very complicated structure when described 
in terms of magnons, will now be configured as the physical, vacuum of Ha 
for a different kind of elementary particle called the spinonH3 

A useful way to characterize the new physical vacuum is through the 
perfectly regular array (|2l]) of Bcthc quantum numbers as illustrated in the 
first row of Fig. H. The spectrum of Ha can then be generated systemati- 



-N/4+1/2 U N/4-1/2 



noioioaoioaoaoi 



Figure 4. Configurations of 
Bethe quantum numbers Ii for 
the N = 32 ground state and 
for one representative of the 
2-spinon (triplet) excitations 
with Si, = St = 1 ■ Each gap in 
o»o«o»o«o»o o • o • one ^-configuration (big circles) 
represents a spinon. 

cally in terms of the fundamental excitations as characterized by elementary 
modifications of this vacuum array. 

In the subspace with = 1, a two-parameter set of states is obtained 
by removing one magnon from the state \A). In doing so we eliminate one 
of the N/2 Bethe quantum numbers from the set in the first row of Fig. [| 
and rearrange the remaining Ii in all configurations over the expanded range 
\h\ < \N. Changing 5|, by one means that the Zj switch from half-integers to 
integers or vice versa. The number of distinct configurations with -fj+i— Ii > 1 
is N(N+2)/8. A generic configuration consists of three clusters with two gaps 
between them as shown in the second row of Fig. |[ The position of the gaps 
between the /j-clusters determine the momenta k\,k2 of the two spinons, 
which, in turn, add up to the wave number of the two-spinon state relative to 
the wave number of the vacuum: Q = k—kA = k\ + k2. 

A plot of the 2-spinon energies E — Ea versus wave number k — kA for 
N = 16 as inferred via ( [l9| ) from the solutions {zi} is shown in Fig. |^ (•). The 
dots, which represent the corresponding data for N — 256, produce a sort of 
density plot for the 2-spinon continuum which emerges in the limit N—r-^ oo. 
The exact lower and upper boundaries of the 2-spinon continuum areEJ 

= |j|sinQ|, uu{Q) =7rJ|sin||. (23) 

Like magnons, spinons carry a spin in addition to energy and momentum. 
Unlike magnons, which have spin 1, spinons are spin-1/2 particles. For even 
N, where all eigenstates have integer- valued S^, spinons occur only in pairs. 
The spins si, S2 = ±1/2 of the two spinons in a 2-spinon eigenstate of Ha can 
be combined in four different ways to form a triplet state or a singlet state. 
We have already analyzed the spectrum of the 2-spinon triplet states 
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Figure 5. 2-spinon 
(triplet) excitations 



with 55 



j t — St — r 
for N = 16 (•) and 

JV = 256 (dots). 
Also shown are the 
2-spinon (singlet) 
excitations with 



2V 



= S T 
16 (o). 



for 



with si = S2 = +1/2 (St = 1, S T = +1). The 2-spinon triplets with 
si = S2 = —1/2 (St = 1, S T = —1) are obtained from these states by a spin- 
flip transformation. One set of 2-spinon states with si = —S2 are the triplets 
with S T = 0. They are obtained via Bethe ansatz by a simple modification of 
the {Ii} configurations from the kind shown in the second row of Fig. || and 
yields an equal number of (real) solutions. Symmetry requires that all three 
triplet components (S T = 0, ±1) are degenerate. 

The 2-spinon singlet states (S T — St — 0) are characterized by one 
pair of complex conjugate solutions z\ = z£ in addition to the real solutions 
Z3, . . . , z N / 2 - Finding the /^-configurations for a particular set of eigenstates 
with complex solutions and then, solving the associated Bethe ansatz equa- 
tions is, in general, delicate task-uiHld The 2-spinon singlets for N = 16 and 
< Q < 7r are shown as open circles in Fig. [H| As N grows large, the effect 
of the complex solutions z\ — z\* on the energy relative to that of the real 
solutions Z3, . . . , zn/2 diminishes. In the limit N — > 00, the 2-spinon singlets 
also form a continuum with boundaries (|2^). Yet the effect of the complex 
solutions will remain strong for other quantities including selection rules and 
transition rates. 

The 2-spinon triplets play an important role in the low-temperature 
spin dynamics of quasi-lD antiferromagnetic compounds such as KC11F3, 
Cu(C 6 D 5 COO) 2 ■ 3D 2 0, Cs 2 CuCl 4 , and Cu^iH^^O^). They are 
the elementary excitations of Ha which can be directly probed via inelastic 
neutron scattering rB The 2-spinon singlets, in contrast, cannot be excited 
directly from \A) by neutrons because of selection rules. The singlet excita- 
tions are important nevertheless, but in a different context. Some quasi- ID 
antiferromagnetic compounds like CuGeO^ are susceptible to a spin-Peierls 
transition, which involves a lattice distortion accompanied by an exchange 
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dimerization.Ha The operator which probes the dimer fluctuations in the 
ground state of Ha couples primarily to the 2-spinon singlets and not at 
all to the 2-spinon triplets. 

5 Spinons Excited from Antiferromagnetic State 

Although the computational application of the Bcthc ansatz yields the exact 
wave functions of the 2-spinon states for very large systems with little effort, 
the very structure of the Bethe wave function makes it hard to use this knowl- 
edge for the calculation of the transition rates pertaining to any fluctuation 
operator of interest. Nevertheless, there exist ways to extract usefuliineshape 
information from the Bethe wave function for specific situations.cH Here we 
discuss an alternative approach, which exploits the higher symmetry of Ha in 
the limit N — ► oo, described by the quantum group Uq^slz)- For example, the 
asymptotic degeneracy of the 2-spinon triplets and singlets noted previously 
is attributable to this higher symmetry. 

The algebraic analysis of a completely integrable spin chain employed for 
the purpose of calculating dynamic structure factors for soecific fluctuation 
operators requires the execution of the following programJla 

• Span the infinite-dimensional physically relevant Hilbert space in the form 
of a separable Fock space of multiple spinon excitations. 

• Generate the eigenvectors of the Hamiltonian in this Fock space by 
products of spinon creation operators (so-called vertex operators) from the 
ground state (physical vacuum). 

• Determine the spectral properties (energy, momentum) of the multiple- 
spinon states accessible from the ground state by the selected fluctuation 
operator. 

• Express the fluctuation operator of interest in terms of vertex operators. 

• Evaluate matrix elements of products of vertex operators as are needed 
for the selected fluctuation operator in the spinon eigenbasis. 

In the following, we outline how this program was carried out for the 
calculation of the exact 2-spinon part of one dynamic structure factor for the 
ID s = \ XXZ model, 

OO 

HxXZ = -J ^2 ( S n S n+l + S n S n+l + AS n S n+l)- ( 24 ) 
n— — oo 

At T — Hxxz has a ferromagnetic phase for A > 1, a critical phase (spin- 
fluid) for — 1 < A < 1, and an antiferromagnetic phase for A < —1. The 
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algebraic analysis operates in the massive phase stabilized by Neel LRO at 
A < — 1, but the isotropic limit A— > — 1 _ , which is equivalent to Ha, can be 
performed meaningfully. 

Each spinon in the m-spinon eigenstate |£ m , e m ; . . . ; £1, ei)j is character- 
ized by a (complex) spectral parameter of unit length and a spin orientation 
e;=±l. The spectral properties then follow from the rules-governing the ap- 
plication of the translation operator and the Hamiltonianiia 

T\^ m ,£m,- ■ = I I 77 > |gm,e m ;- ■ -^l^l) 1-J : (25) 

m 

H\€m,erri- ■ J = ^ e (&)l£m,eml- ■ -^l^^i, (26) 



r(0=e -^) =r xi^ ) e(0 = J^4logr(0, (27) 



6 x (y) = (x;x)(y;x)(xy- 1 ;x), (y; x) = JJ (1 - yx n ). (28) 

Here 5 is the deformation parameter of the quantum group U q (sl2)- Its value 
is determined by the exchange anisotropy, A = (q + q~ 1 )/2. The twofold 
degenerate vacuum is represented by the vectors |0)o,|0)i, which transform 

into each other via translation: T\0)j = 1 0) 1 j . In the Ising limit A — > —00, 

they become the pure Neel states | . . . fj.fl ■■■);!•■■ itlf • • ■ )• 

By eliminating £ from the relations (|2^) and ( |26| ) we obtain the spinon 
energy-momentum relation 

e x (p) = - k 2 cos 2 p, 1 = sinh — — , (29) 

which is independent of the spin orientation. The elliptic integrals K = 
K(k), K' = K(k') and their moduli k, k' = yl — k 2 depend on the anisotropy 
via — q = exp(—TrK'/K). 

The of 2-spinon spectrum with energies #(£1,62) = e (£i) + e (&2) and 
momenta P (£,1,^2) = p(£i) +^£2) is then a two-parameter set with fourfold 
spin degeneracy £1,62 = ±1.E3 In the isotropic limit, they are the 2-spinon 
triplets and singlets discussed in Sec. [|. In a finite system, the singlet-triplet 
degeneracy is removed (see Fig. ^|). Exchange anisotropy splits up the triplet 
levels as well. 
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Q/K 



(%-QJAl 



Figure 6. 2-spinon exci- 
tation spectrum for k = 
0.99 (A ~ -2.305). It 
consists of two partly 
overlapping sheets C_ 
and C+. Sheet C+ 
lies between loo(Q) and 
uj— (Q) in the range 
Q K < Q < n - Q K 
and between u>+(Q) and 
oj-(Q) for 7r - Q K < 
Q < n. The in- 
terval (0, 7r) represents 
one half of the extended 
Brillouin zone. 



Q/n 



In the (Q,lu) plane, each set of 2-spinon excitations forms a continuum 
(see Fig. ||) of two sheets C± with boundaries 

MQ) = -^-sinQ, w±(Q) = + k 2 ± 2kcosQ, (30) 

1 + ft 1 + ft 

where k = cosQ K = (1 — fc')/(l + fc') is a convenient anisotropy pa- 
rameter. The evaluation of the dynamic structure factor of any fluctua- 
tion operator Fq, Eq. ([|), requires that we know exact matrix elements 
j(0|A;|£ m , e m , . . . ,£i! e i)j' for the associated local operators. In principle, 
they can be calculated exactly if we are able to express Ai in terms of vertex 
operators. In practice, the exact results available thus far are limited to the 
dynamic spin structure factors 

+ OQ poo 

Sw(Q,u)= £ / dte^+^iSZWSS), fi = x,z. (31) 

In the m-spinon expansion S ZZ {Q, uj) — S Z °J(Q , uj) + si~JQ , uj) + . . . we 
know the matrix elements that determine the leading termJ13 

M0\S z n \0)j = ( ^^) 2 (- 1 )" +J = 2m„ (32) 

yielding ^(ftw) = Ah 2 5{uj)5(Q — 7T)m z , where m z is the staggered 
magnetization.EJ No exact transition rates have been evaluated for the 2- 

(21 

spinon part S zz (Q,uj) or any higher order term in the m-spinon expansion. 

The corresponding m-spinon expansion of S xx (Q,uj) starts with m = 2, 
because the operator S% does not connect the vacuum sector with itself. All 
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non- vanishing matrix elements which are needed for S x J(Q,uj) can be ex- 



pressed, as it turns out, by a single function JL J (£2 , £1) = j( 
which was determined by Jimbo and Miwa:E£l 



2 (g 4 ;g 4 ) 2 (-gte2) 1 - j ZMe 1 /8)QA-ZT 2 & 2 'f j ) 



- (<?w) 3 M-erV)M-sV) 

where (x; y; z) = n^ m =o ( x _ xy n z m ) and 

(<? 4 £;<? 4 ;<7 4 )(rV;<? 4 ) 



<4 16,-^1, 

(33) 



7(0 



^ ;<7 4. 9 4 )((?2r i ;g 4 ;g 4 



(g 4 ;g 4 ;g 4 ) 
(<? 6 ;<zW)' 



These ingredients yield the following expression to be evaluated:^ 



x{6[Q + P{Si,&)]\x fa,&)+xHt2,€i)\ 2 

+S[Q -tt + Pfa,t2)]\X°fa,ti) - X'&^i)] 2 }- 
A compact rendition of the exact result reads0 

/ - /,-^,- fk.{P-) I tan(Q/2)|- 



c=± 



where 



U 2 



± cosQ 



T = 



^(Q,c) = l±^F 



arcsm 



2IujW c 
k(1 + k)wq 



1^09) =exp 



1=1 



e 7 ' cosh(27/) cos(tfl) - 1 
sinh(2/7) cosh(7Z) 



(34) 



(35) 



(36) 



(37) 



(38) 



(39) 



7 = ttK'/K, t = 2/3/K', and fldix) is a Neville theta function. This function 
is plotted in Fig. [7] for (Q,lu) e C + . S x 2 J(Q,uj) has a square-root divergence 

at the portion ojq(Q) of the spectral threshold. Along the portion uj + (Q) of 

(2) 

the lower boundary and along the entire upper boundary of C+, S X x{Q,(jj) 
has a square-root cusp. The line shapes for (Q,uj) € C— are inferred from the 



symmetry property SxJ(ir — Q,w) = S xx (Q,u) 



(2), 
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The isotropic limit A — » — 1~ is delicate because of its singular nature. 
As the LRO in the spinon vacuum vanishes, the size of the Brillouin zone 
changes from (— 7r/2, +tt/2) to (— w, +ir). A practical consequence of this 
phase transition is that we switch our perspective from considering both sheets 
C± of 2-spinon excitations over the range (— n/2, +7r/2) to considering only 
the sheet C+, now with boundaries (p3|), over the extended punge (— ir, +ir). 



With these subtleties taken into account, Eq. (36) reduces tc£3 



5( 2 j(g,^) = i[^(Q)-c 2 ]- 1 /2 e -/W j (40) 



OC 



cosh 2x cos xt — 1 x irt u _ t j Ujfj(Q) — lo 2 l (Q) 



m= h ^ *sinh2* CO shx ^' T^ C ° Sh V »>-u>l(Q) • ^ 

This function, which is plotted in Fig. ^, has a square- root cusp at uju(Q) and 
a square-root divergence with logarithmic corrections at ujl{Q)- 
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Figure 8. Exact 2- 
spinon dynamic spin 
structure factor. The 



nonzero 



expression 
only in the shaded re- 
gion of the (Q.oj)-plane 
bounded by lol(Q) and 
£Jr/(Q). 



6 Solitons Excited from Neel State 



Near the Ising limit (A — > — oo), the exact result (36) can be expanded in 
powers of the anisotropy parameter k. To leading order, we obtain 



1 



2 cos 2 Q 



/cos 2 Q- ( - ) (-)(/,-|c. ^()\ - o|) 



(42) 



with = uj/2I — 1, which is identical (in that order) to the result obtained 
by Ishimura and ShibatJ from a first-order perturbation calculation about 
the Ising limit. In their calculation, the two vacuum states are approximated 
by thejjure Neel states (see Fig. ||). The spinons are kink solitons (domain 
walls) ,0 which produce two adjacent up or down spins in an otherwise un- 
perturbed Neel configuration. The only states which contribute to S xx {Q,uj) 
in leading order of the perturbation calculation are states which contain two 
solitons with equal spin orientation. The states with solitons of opposite spin 
orientation dominate the spectrum of S zz (Q,uj) except for the central peak 
at u) = 0, Q = 7r, which involves a transition between the two vacuum states. 



Figure 9. Local basis 
states from which trans- 
lationally invariant Neel 
states and 2-soliton states 
are constructed. The 
latter are limiting cases 
of the 2-spinon states 
discussed in Sec. ^. 

The 2-soliton spectrum with boundaries < kcosQ, is akin to the 



vacuum (Sr =0) 
2-soliton (S?=l) 

2-soliton (S^=0) 



■j> imtmtmt) 

i i> uti0iti(Q)it> 

l®Jtlt(0)tl> ut(0)tlll(Q)lt> 



■ t> 
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generic 2-spinon spectrum shown in Fig. |6| but it does not capture any of the 
features that characterize the regime Q K < Q < n — Q K around the Brillouin 
zone boundary. The limitations of the perturbation results are much less 
severe near the center of the Brillouin zone. This illustrated in Fig. [Hi which 
compares the line shapes of Eq. ( |36| ) and (^) at Q — n. Here a fairly 
strong departure from the Ising limit is required before a significant difference 
between the two results is produced. 



Figure 10. Line 

shapes of S^J(Q,co) 
at the zone center 
Q = it as established 
by the exact result 
(Bq) (solid line) and 
predicted by the first 
order jwturbation 
result (||) (dashed 
line). 



0.5 1 1.5 2 

CO// 

Two intensively studied physical realizations of Hxxz at A < —1 arc 
CSC0CI3 and CsCoBi-3. Spectroscopic data which probe both the quantum 
and thermal fluctuations-of those materials are available from several neutron 
scattering experimentsEa 
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